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ABSTRACT 

The X-ray emission from active galactic nuclei (AGN) is highly variable. Measurements 
of time lags (characterised by lag spectra) between variability in the light curves in 
energy bands corresponding to directly observed continuum emission from the corona 
around the black hole and to X-rays reflected from the accretion disc adds a fur- 
ther dimension to studies of the structure and energetics of these systems. We seek 
to understand these measurements in terms of the physical parameters of the X-ray 
source (its location, extent, etc.) through the calculation of theoretical lag spectra for 
a range of source parameters in general relativistic ray tracing simulations, combined 
with knowledge of the observed variability of the X-ray emission from AGN. Due to 
the proximity of the emission to the central black hole, Shapiro delays are important 
and the effects of general relativity should be considered when interpreting the lags as 
the light travel time between the source and reflector. We show that it is important 
to consider dilution of the lag by the contribution of both the primary and reflected 
spectral components to the observed energy bands rather than observing pure contin- 
uum and reflected emission, reducing the measured lag by up to 75 per cent compared 
to the 'intrinsic' time lag due to light travel times. We find that the observed lag 
spectrum of the narrow line Seyfert 1 galaxy 1H 0707-495 implies an X-ray source 
extending radially outwards to around 35 r g and at a height of around 2 r g above the 
plane of the accretion disc, consistent with the constraints obtained independently by 
considering the emissivity profile of the accretion disc. By investigating the influence 
of the propagation of X-ray luminosity fluctuations through the source region we find 
it is possible to reproduce the shape of the low frequency part of the lag spectrum 
(where the hard 'primary' band lags behind the soft 'reflected' band) as the effect of 
luminosity fluctuations originating in the centre of the X-ray source, close to the black 
hole, and propagating outwards. 

Key words: accretion, accretion discs - black hole physics - galaxies: active - X-rays: 
galaxies. 



1 INTRODUCTION 

High resolution X-ray spectra from long observations of ac- 
tive galactic nuclei (AGN) have yielded an unprecedented 
amount of information about the structure of these systems, 
detecting emission from material right down to the inner- 
most stable orbit and even the event horizon probing not 
only the mechanisms at work in AGN but also the strong 
gravity regime close to the black hole itself. 

The observed X-ray spectra can be analysed in the con- 
text of emission from a primary X-ray source in a corona 

* E-mail: drw@ast.cam.ac.uk 



of energetic particles close to the central black hole, which 
inverse-Compton scatter the thermal photons emitted from 
the accretion disc up to X-ray energies. T his emission is fol- 
lowed by reflection from an accretion disc (|George fc Fabianl 
ll99lT) which is g eometrically thin in the eq uatorial plane and 
optically thick (IShakura k, Sunvaevll973h . Detailed analysis 
of reflection features in the spectrum, notably the Kq emis- 
sion line of iron at 6.4 keV, has allowed constraints to be 
placed on the properties inc luding the location and exten t 
of the primary X-ray source |Wilkins fc Fabi an 2011, 20121 '). 

The X-ray emission from AGN is highly variable and 
recent studies of this variability have added a further di- 
mension to the studies of these systems. The narrow line 
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Figure 1. The lag spectru m from a 1 Ms lightcurve of 1H 0707- 
495 from lKara et aL showing the time lag for the variabil- 

ity at different Fourier frequencies between the hard (1.0-4.0 keV) 
band corresponding to the primary continuum emission and the 
soft (0.3-1.0 keV) band corresponding to reflection from the ac- 
cretion disc. A positive lag indicates that variability in the soft 
(reflected) band lags behind that in the hard (primary) band. 



Seyfert 1 galaxy 1H 0707-4 95 was the first AG N found to 
exhibit reverberation lags (jFabian et alj [20091 ) . where the 
variability in the spectral band dominated by reflection from 
the accretion disc is found to lag behind that in the spec- 
tral band corresponding to direct emission fro m the coronal 
X-ray source by tens to hundreds of seconds (|Zoghbi et al.l 
2010). Since then, reverbera tion lags have been discovered 
in a multitude of othe r AGN lEmmanoulopoulos et al.ll201ll; 
de Marco et al.ll201ll; IZoghbi fc Fabian! l201ll ; IZoghbi et"al] 
2012l : lde Marco et al.ll2012| ) and even in Galactic stellar mass 
X-ray binaries where the black hole mass and length scale 
invariance gi ves equivalent reve rberation lags of the order 
milliseconds (|llttlev et al.l[201lh . 

Measuring these phase lags measures the light travel 
time between the source and reflector and allows us to probe 
scales as small as 10 light-seconds in AGN. The reverbera- 
tion of X-rays fr om the accretion d isc is characterised by 
the lag spectrum (|Nowak et al.|[l999l ) which shows the time 
lag of variability at different frequencies between the energy 
bands. The lag spectrum o f 1H 0707-495 was fi rst obtained 
from a 500 ks lightcurve bv lZoghbi et~ai1 (|2010f ) and then in 
more detail followi ng observations with XMM Newton to- 
talling 1.3 Ms by iKara et "ail <|2012r i. Such a lag spectrum, 
from a 1 Ms lag spectrum of 1H 0707-495 between the 'hard' 
(1.0-4.0 keV) band dominated by the primary continuum 
and the 'soft' (0.3-1.0 keV) band dominated by reflection 
from the accretion disc is shown in Fig. [1] 

At temporal frequencies of around 10 -3 Hz, a lag of 
around 30 s can be seen between the arrival of the X-ray con- 
tinuum and the reflection from the accretion disc. Naively, 
this can be converted into the distance of the reflecting 
material from the primary X-ray source, which taking the 
mass of the black ho le in 1H 0707-495 to be 2 x 10 6 M Q , 
i|Zhou fc Wand 120051 ), gives a distance of 2r g (a gravita- 
tional radius, lr g = ^-^-). This simple calculation, however, 
takes into account neither the multiple paths that X-rays 
may follow from the source to the reflector nor the effects 
of general relativity which are expected to be significant in 
such close proximity to the black hole. It is assumed that the 



spectral energy bands correspond to detections of pure con- 
tinuum and reflected emission, while in reality there will be 
contributions from both spectral components in each energy 
band. Finally, this analysis does not explain the full shape 
of the lag spectrum, rather it just looks at the longest time 
lag where the reflection follows the primary continuum. 

The reverberation of of the iron Ka emission line from 
the accretion disc of an AGN due to a localised flare of 
X-ray emission fro m the c orona was fi rst considered by 
i|Fabian et al.lll989T ) and by IStellal (Il990l ) in the context of 
measu ring the mass of the central black hole. iGilfanov et alj 
(2000) quantified the effect of reverberation in the Galactic 
black hole binary Cygnus-Xl, computing an approximate 
solutio n neglecting relativistic effects while iRevnolds et al.l 
( 1999) accounted for the transport of the observed X-rays 
in general relativity to compute '2D transfer functions' giv- 
ing the observed flux as a function of observed photon en- 
ergy as time elaps e s from a single, localised, X-ray flare. 
lYoimg fc Reynolds! l|2000T) related this work to observable 
signatures that may be detected by future generation X-ray 
missions in the context of comparing observed reverberation 
of X-ray flares with a library of computed transfer functions 
to determine the mass of the black hole as well as the loca- 
tion of the X-ray flare. 

We here present an analysis of X-ray reverberation ap- 
plicable to the current state-of-the-art observations with 
XMM Newton. A systematic analysis of lag spectra com- 
puted theoretically for reflection of X-rays originating from a 
variety of X-ray sources of different sizes and geometries and 
in different locations allows us to understand the observed 
form of the lag spectra in AGN and to ask how much can 
be learned from these measurements about the energetics 
and geometry of these systems. While applicable generally 
to X-ray sources exhibiting a reverberation lag, this study 
is motivated by the observed lag spectrum of 1H 0707-495. 
As well as the location and geometry of the X-ray source, 
we consider observational constraints in separating the de- 
tection of the primary continuum and the reflected X-rays 
from the accretion disc and we consider how variations in 
luminosity may propagate through the X-ray source region. 



2 LAG SPECTRA 

Following iNowak et alj (|l999l ). time lags between spectral 
components (with their own light curves) are characterised 
by the lag spectrum, showing the phase (and hence) time 
lag between the Fourier frequency components that make 
up the light curves. The light curve, F(t), is considered to 
be the sum over components of all frequencies: 



F(t) = / F(uj)e iuJt du 



The amplitude and phase of each frequency component 
is given by the Fourier transform of the light curve and can 
be written F(w) = -F(w) e %v and is computed by 

F{u) = f F{t)e- l ^dt 



The phase lag between two spectral components, say the 
hard and soft light curves, H(t) and S(t) respectively, can 
be found by considering the complex form of their Fourier 
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Figure 2. Lag spectrum between a primary light curve (taken 
from the light curve of 1H 0707-495 in the 1-4 keV band) and the 
same light curve delayed by 30 s. The lag spectrum shows a con- 
stant lag of approximately the input 30 s across all frequency com- 
ponents until phase- wrapping occurs when the lag corresponds to 
a half- wave shift in phase at / = i . The drift to slightly a longer 
lag at the lowest frequencies is due to the error in the Fourier 
transform of the observed light curve. 

transforms H(ui) = |fl'(w)|e <v ' and S(u) = e if? and 

computing the cross spectrum 

G(u>) = S*(oj)H(lu) = |5H| \h(uj)\ e liv ~ e) (1) 

The argument of which gives the time lag, r, since tp = cut 
(and converting from angular to linear frequency, /) 

r(/) = ^arg(c(/)) (2) 

Following this calculation (and sign convention), a positive 
time lag indicates that the variability in the soft energy 
band, S(t), is lagging behind that in the hard energy band, 
H(t), as would be expected if the hard band is dominated by 
an X-ray continuum which is reflecting off of the accretion 
disc and this reflected radiation dominates the soft band. 

We here issue a warning to the reader that this defi- 
nition is the oppos ite s ense of the time lag as adopted by 
IZoghbi et all l|20ld ) and lNowak et all l| 19991 ) who define the 
soft band lagging behind the hard band to be the negative 
lag, though defining this as a positive lag is more instructive 
when considering the features of lag spectra arising from re- 
verberation of X-ray variability from a reflecting accretion 
disc. 

Detection of a time lag between spectral components re- 
quires the light curves to be coherent (i.e. there is a constant 
phase lag between the Fourier components over a given range 
in the two light curves, rather than a random one), although 
iKara et al.l |2012l ) argue that total coherence is not required 
as uncorrelated (incoherent) variability will not contribute 
to the calculated lag as this part will average to zero. 

To explore the general features of reverberation lags, the 
lag spectrum was computed for the simplest case of a 're- 
flected' light curve exactly following the primary light curve 
but delayed in arrival by 30s (i.e. shifting the time axis by 
30 s). Such a lag spectrum, using the light curve of 1H0707- 
495 in the 1-4 keV band as the input is shown in Fig. [2] 

The lag spectrum shows a constant lag of approximately 
the input 30 s across all frequency components until the lag 
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corresponds to a half-wave shift in phase at / = -=- (give or 
take the error introduced by the numerical, discrete Fourier 
transform of the real light curve and the temporal binning 
of the light curve, the combination of which causes the drift 
to a slightly longer lag at the lowest frequencies). Once / = 
the waveform could have been shifted either forwards 
or backwards by half a wave and since the phase angle from 
which the lag is calculated is defined to be in the range — n < 
tp < n, phase wrapping occurs. At higher frequencies, the 
phase continually wrap causing the lag to oscillate around 
zero until it decays away. It is not possible to measure lags 
longer than a shift of half a wavelength. 



3 SIMULATING LAG SPECTRA FOR X-RAY 
REFLECTION 

3.1 The Transfer Function 

The time lag between the observation of a change in the 
primary X-ray continuum and the corresponding change in 
the reflected flux from the accretion disc is characterised 
by the transfer function, T(t). The transfer function is the 
response seen from the accretion disc to an instantaneous 
flash of light from the primary X-ray source. 

The transfer function is calculated by general relativis- 
tic ray tracing simulations in two stages. Initially, rays are 
transferred from the source to the accretion disc, in the equa- 
torial plane of the Kerr geometry around the central black 
hole with dimensionless spin parameter, a, as described in 
IWilkins fc Fabian! |2012T l. 

Initially, the X-ray source is assumed to be an isotropic 
point source, emitting equal power into equal solid angle 
in its own instantaneous rest frame. The photon's initial 4- 
momentum is transformed from the source's frame, a locally 
flat spacetime described by a tetrad of orthonormal basis 
vectors, to the global Boyer-Lindquist co-ordinate system to 
calculate the constants of motion of the ray, k, h and Q. 
The ray is then traced by numerical integration of the null 
geodesic equations (the equations of motion of a photon) in 
the Kerr spacetime. Taking c = 1 and = 1 to work in 
natrual units of r g = the geodesic equations can be 

written: 

[(r 2 + a 2 cos 2 0)(r 2 + a 2 ) + 2a 2 r sin 2 0] k - 2arh 

r 2 (l + " 2 c r T 2 6 - |) (r 2 + a 2 ) + 2a?r sin 2 

2ark sin 2 + (r 2 + a 2 cos 2 -2r)h 
* ~ (r 2 + a 2 )(r 2 + a 2 cos 2 - 2r) sin 2 + 2a 2 r sin 4 
■ Q + (ka cos — h cot 0) (ka cos + h cot 6*) 

7 

r 2 = — \ki — hip — p 2 # 2 J 

Where the dots denote differentiation with respect to some 
affine parameter, a. 

Spatially extended X-ray sources may be explored 
through Monte Carlo ray tracing simulations. Rays are 
started at random locations (with a uniform probability 
density function such that all locations are equally likely) 
within an allowed cylindrical source region (defined by a 
lower and upper height from the disc plane as well as an in- 
ner and outer radius) and are assigned random initial direc- 
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tion cosines cos a and /?, again with uniform probability den- 
sity. This will simulate the effect of an X-ray source of finite 
spatial extent that is optically thin to the X-rays it emits. 
Each local region of the source is taken to be co-rotating 
with the element of the accretion disc in a relativistic Kep- 
lerian orbit at the same radius. This model does, however, 
assume the surface brightness over the entire extent of the 
source varies simultaneously, which may be unphysical if the 
source region extends across multiple gravitational radii, as 
the light crossing time of the source is no longer negligible 
compared to the light travel time from the source to the 
reflector, however this model may be used to approximate 
the behaviour of an extended source and we will revisit the 
propagation of fluctuations in SJU 

When the rays reach the accretion disc, their arrival 
time as measured by an observer at infinity (the t co- 
ordinate) as well as position (r, cp) are recorded and they 
are sorted into radial and azimuthal bins on the accretion 
disc. 

The travel time of photons from each radial and az- 
imuthal bin on the accretion disc as well as their energy 
measured by an observer at infinity is determined by tracing 
parallel rays, travelling perpendicular to a flat image plane 
a large distance from the black hole, representing the area 
of the sky imaged by the X-ray telescope. Rays are started 
on a regular grid (the 'image plane') travelling perpendicu- 
lar to this plane plane, a distance 10000 r s from the black 
hole (to be certain the spacetime here is flat), inclined to 
the rotation axis of the black hole at the angle at which the 
telescope is observing the system, henceforth taken to be 
53deg, the measured inclination of 1H 0707-495. The image 
plane is centred on the black hole and spans the area of in- 
terest (i.e. the accretion disc, assumed to have a radius of 
400 r g ). 

Due to the symmetry of the Kerr metric, propagating 
rays backwards in time (starting at the image plane to work 
out where they originated from on the accretion disc) is 
equivalent to propagating forwards in time, having reversed 
the direction of spin of the black hole. Propagating the rays 
from the image plane with the spin of the black hole reversed 
until they reach the equatorial plane will therefore give the 
position on the accretion disc (r and ip co-ordinates) from 
which the ray originated as well as the travel time of the 
ray from the disc to the telescope as measured by the ob- 
server at infinity (t co-ordinate). The rays are sorted into 
the radial and azimuthal bins upon the accretion disc and 
the average travel time of a ray is calculated for each bin. 
The number of rays landing in each bin is counted, which 
gives the area of that bin projected onto the image plane by 
general relativistic light bending and will determine the flux 
received from each part of the accretion disc. 

The energy and arrival rate of the photons arriving 
along a given ray measured by the observer at infinity is 
computed by the projection of the photon 4-momentum to 
the observers' timelike axes (which are their 4- velocities), 
though before doing so, since the photon under considera- 
tion is propagating backwards, the spacelike components of 
the 4-momentum are reversed (so that the photon is trav- 
elling in the correct direction with respect to the orbiting 
material in the accretion disc and the correct Doppler shift 



is computed). 

-l _ "a = vq ■ p(Q) = g^oP u {Q) 

~ W ve-p(E) g pa v p E p°(E) U 

The average energy shift for each bin on the accretion disc 
is computed from all the rays originating from it. It is as- 
sumed that the primary X-ray source is emitting photons of 
all energies and that this will excite emission lines and other 
spectral features from the disc. The energy shift of the ob- 
served photons is assumed to arise only from the transport 
of rays from the accretion disc to the observer, taking into 
account special relativistic Doppler shifts as well as gravita- 
tional redshifts. Material in the accretion disc is assumed to 
be in a (relativistic) prograde Keplerian orbit in the equa- 
torial plane around the black hole. 

When a photon from the primary X-ray source lands in 
a specific bin on the accretion disc, the time from the disc to 
the observer as well as the energy of the photons measured 
at infinity is looked up from the values calculated for that 
bin from the rays traced from the image plane to the disc. 
The total arrival time of the ray from the primary X-ray 
source, via the disc, to the observer is the sum of the two 
computed times. 

Finally, the transfer function is computed by counting 
the number of rays that arrive at the observer as a function 
of time, N(t, r, ip, E), integrating over all energies and posi- 
tions on the disc. Assuming the X-ray source emits photons 
at an equal rate along each ray in its own rest frame (i.e. it 
is isotropic), the arrival rate of photons along each ray mea- 
sured by the observer will vary according to the redshift, 
^ = g , which is the product of the redshifts from the 
source to the disc (g^) and from the disc to the observer 

tec,)- 

T(t) = j N(t, r, p, E) gsngoh rdrdipdE (4) 

The ray tracing code to compute these transfer func- 
tions was implemented to run rapidly on graphics processing 
units (GPUs) using the nvidia CUDA framework, allowing 
the paths of hundreds of rays to be computed in parallel (for 
more details of the impleme ntation of general relativ istic ray 
tracing codes on GPUs, see IWilkins fc Fabianll2012l ). 

3.2 The Arrival Time of the Primary Continuum 

The ray tracing simulations described above give the pho- 
ton arrival rate as a function of time since the initial flash of 
light left the primary X-ray source. However, when measur- 
ing reverberation in observations of AGN, it is the time lag 
between the arrival of the primary continuum and the re- 
flected rays that is measured. One must therefore know the 
arrival time of the primary X-ray emission at the observer 
to define the zero in the time series of the transfer function. 

For an infinitesimal point source, exactly one ray emit- 
ted from the source will end up at any given observer at 
infinity. Therefore the arrival time of the primary X-ray con- 
tinuum may be computed by tracing one ray from the source 
to the observer. Table[T]shows the direct ray travel times (as 
measured by an observer at infinity) from point sources at 
varying heights above the black hole to an observer at a dis- 
tance of 10000 r g from the black hole inclined at 53 deg to 
the rotation axis. A source above the plane of the accretion 
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Table 1. Travel time (in natural units) as measured by an ob- 
server at infinity for a ray propagating from point sources at vary- 
ing heights above the black hole on the rotation axis to an ob- 
server 10000 r g away inclined at an angle of 53 deg to the rotation 
axis. The Shapiro delay is characterised by the difference between 
the travel time of the ray in general relativity and what it would 
be in classical, Euclidean space. 



Source Height 


Ray Travel Time 


Delay 


1.235 r g 


10029 GM/c 3 


30 GM/c 3 


2r g 


10020 GM/c 3 


21 GM/c 3 


5r g 


10014 GM/c 3 


17 GM/c 3 


10 r g 


10008 GM/c 3 


14 GM/c 3 


20 r g 


10001 GM/c 3 


13 GM/c 3 



disc should be closer to the observer than the black hole and 
classically, in natural units, light takes l^p- to travel lr g . 
These direct rays taking longer than 10000 to reach the 
observer illustrates the Shapiro delay here where the pas- 
sage of light is slowed down a s it travels in the proximity of 
the black hole (|Shapirolll964r ). 



3.3 The Reflected Light Curve and Lag Spectrum 

The transfer function describes the response seen in the re- 
flection from the accretion disc due to an instantaneous flash 
of X-rays from the primary source. Hence, if the primary 
X-ray source, observed in a 'hard' band for the X-ray spec- 
trum is varying as described by some light curve H(t), the 
response in the reflection seen in a 'soft' band, S(t), is given 
by the convolution of the primary light curve with the trans- 
fer function. 

S(t) = H(t)®T{t) = J H(t')T(t-t')dt' (5) 

Ray tracing in the simulations is carried out in natural units, 
measuring time in units of . The time series of the trans- 
fer function is trivially converted to seconds for the mass of 
the black hole of interest. 

In order to generate a realistic lag spectrum with the 
power spectral density (PSD) observed in the primary con- 
tinuum of an AGN, the primary, hard light curve is taken 
to be the observed light curve of the narrow line Seyfert 1 
galaxy, 1H 0707-495 in the 1-4 keV band in 10 s bins. Using 
the observed light curve of 1H 0707-495 also means that lag 
spectra are simulated with a realistic time resolution that is 
achievable with the current generation of X-ray telescopes. 
The mass of the black hole is t aken to be 2 x 10 6 M q , as 
quoted in the literature (see e.g. IZhou fc Wang |2005). 

The primary light curve is convolved with the trans- 
fer function obtained from the ray tracing simulations in 
order to obtain the light curve of the reflected component. 
The Fourier transforms of both the primary and reflected 
light curves were then computed. The cross and lag spectra 
are calculated using Equations [T] and [5] The convolution, 
Fourier transforms (using the Fast Fourier Transform, FFT 
algorithm) and lag calculations were performed in matlab. 



3.4 Extended Sources 

For an extended source, a range of rays following different 
paths will be observed by the telescope as the primary con- 
tinuum. A transfer function (photon arrival rate vs. time) 
for an extended source region may be calculated in exactly 
the same way as the rays are traced from the observer to 
the accretion disc. If the source is said to extend radially at 
some height above the accretion disc, rays are traced from 
the image plane until they reach this plane and their time is 
measured if they intercept the region of the plane spanning 
the source. 

In order to compute the lag spectrum arising from the 
reflection of X-rays originating from an extended source, the 
observed light curve is taken to be the 'intrinsic' variability 
of the X-ray source. This is permissible since the transfer 
functions do not impart any power on the variability spec- 
trum, they just give a time lag and this will mean the sim- 
ulated variability has a physically realistic power spectrum. 
Simulated light curves that would be observed for the pri- 
mary (hard) and reflected (soft) spectral components, H(t) 
and S(t) respectively, are computed by convolving the light 
curve with the calculated transfer functions for each com- 
ponent. 

H(t) = H ohs (t)®T Bt . c (t) (6) 
S(t) =H ohs (t)®T let (t) (7) 

The lag spectrum between the two spectral components is 
calculated as above. 



4 THEORETICAL LAG SPECTRA 

4.1 Lag Spectra for Point Sources 

Initially, the transfer functions and lag spectra were com- 
puted for reflection from isotropic point sources located at 
varying heights on the rotation axis above the black hole, as 
shown in Fig. [3] 

A striking feature of the transfer functions is the peak 
in X-rays after the initial arrival of the reflection from the 
accretion disc. This corresponds to the 're-emergence' of 
the redshifted, gravitationally lensed (and thus magnified) 
reflection from the far side of the accretion disc, travel- 
ling round the black hole to the observer as was noted by 
iRevnolds et~al] j 19991 ). 

Turning to the resulting lag spectra, a constant lag is 
observed at low frequencies, corresponding to the average 
time lag of the transfer function 

i = - F —— f tT(t) dt (8) 
/ T(t) dtj y > y ' 

Towards higher frequencies, the contribution to the lag 
from the extended tail of the transfer function decreases 
as phase wrapping causes these contributions to the lag to 
average to zero. This decreases the average lag time at higher 
frequencies, until the peak in the transfer function is left to 
dominate, causing the lag spectrum to flatten off before the 
jump down where phase wrapping occurs for all components 
of the lag. 

The trans fer functions roughly fo llow the form of those 
computed by iGilfanov et ail ^OOOT l who calculated the 
transfer functions for point sources in Euclidean spacetime, 
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Figure 3. [(a)] Transfer functions and |(b)| lag spectra for reflection due to isotropic point sources at varying heights, stationary on 
the rotation axis above the black hole, mass 2 X 10 6 Mq . The accretion disc is observed at an inclination angle of 53 deg. A positive 
time lag indicates that variability in the reflected X-rays lag behind variability in the primary emission from the source. 



taking into account no relativistic effects and only for ac- 
cretion discs extending inward to 10 r g . Here, the fully rel- 
ativistic transfer functions are computed for accretion discs 
extending all the way down to the innermost stable orbit for 
a maximally spinning black hole. As in the classical case, the 
transfer functions fall of as t~ 2 at long times, though the rel- 
ativistic effects cause the initial peak in the transfer function 
to be steeper. Extending the accretion disc in to the inner- 
most stable orbit means the flatt ened part visib l e at e arly 
times in the transfer functions of iGilfanov et all (|2000h are 
not seen (consistent with their trend as the inner radius of 
the disc is decreased). 

Considering, for illustration, the case of the point source 
at a height of 5 r g , the mean arrival time of rays in the trans- 
fer function is around 290 s after the arrival of the primary 
continuum from the point source. This time can be com- 
pared with travel times computed for the equivalent source 
and reflector in flat, Euclidean spacetime. For the case of 
a point source at 5r g , the primary continuum takes 168 s 
longer to reach the observer accounting for relativistic ef- 
fects than it would in Euclidean spacetime and computing 
the transfer function using Euclidean ray travel times (but 
keeping the relativistic formulation for the projection of ar- 
eas on the accretion disc to the observer and the emissivity 
profile of the disc), the average arrival time is delayed by 
173 s. Since the lag is measured between the arrival of the 
continuum and the reflection, relativistic effects on the light 
travel time increase the measured lag by around 4 s (the 
'Shapiro delay'). When directly equating the lag to the light 
travel time from the source to the reflector, this means the 
source appears 0.5 r s further from the accretion disc than it 
really is. 



4.2 Lag Spectra for Extended Sources 

Whilst it is instructive to consider the lag spectra arising 
from point sources in general relativity, observations of the 
relativistically broadened iron K line profile in 1H 0707-495 
suggest that the X-ray source is rad i ally e xtended out to 
around 35 r s jWilkins fc Fabianll201ll . I2012T I. It is therefore 



important to understand the form of lag spectra arising from 
extended sources. 

While the following analysis is generally applicable to 
X-ray reverberation from spatially extended sources, we are 
motivated by the observed lag spectrum of 1H 0707-495 and 
illustrate the effects with source geometries similar to that 
suggested by the emissivity profile observed in this object. 
Further more, we consider a source region whose luminosity 
is constant throughout its extent. While in reality the lu- 
minosity is likely to vary across the source, perhaps being 
concentrated in the central regions, we here elect to min- 
imise the number of free parameters such that we may gain 
insight into the physical effects. When assuming a constant 
source luminosity, the modelled source extents will represent 
the bulk of the emitting material. 

The transfer functions and lag spectra for reflection of 
X-rays from source regions of constant height and vertical 
extent but increasing in radial extent outward from the rota- 
tion axis are shown in Fig. [4] The most notable feature of the 
transfer function for extended sources is the extended tail 
to long reverberation times since from an extended X-ray 
source, there are many possible light paths for the primary 
X-rays to reach the accretion disc. 

In addition to the extended tail of the transfer function, 
some reflected rays (those originating from the edge of the X- 
ray source closest to the observer reflecting off the part of the 
disc immediately below) may reach the telescope before the 
mean arrival time of the continuum, meaning the transfer 
function begins before the time t — 0. As for the case of a 
point source, the reflected X-rays are seen to rise after their 
initial arrival owing to the re-emergence of X-rays from the 
far side of the accretion disc, lensed around the black hole, 
though this is less pronounced due to the increased range of 
ray paths now possible with an extended source. 

a constant 



Turning to the lag spectrum in Fig. 4(b) 



lag is once again observed at low frequencies corresponding 
to the mean arrival time of reflected rays over the transfer 
function, though with an extended X-ray source, since more 
of the flux is in the extended tail of the transfer function, 
there is a much more pronounced decay of the lag measured 
in the higher frequency components. For the most extended 
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Figure 4. |(a)| Transfer functions and |(b)| lag spectra for the instantaneous illumination of X-ray sources extended vertically between 
2r g and 2.5 r g above the plane of the accretion disc and with different radial extents outward from the rotation axis. 



sources, there is almost a continuous decay of the measured 
lag to zero in the high frequency components with no phase- 
wrapping observed. The extremes of the source being further 
from the black hole reduces the shortest possible light travel 
time to the observer since further from the black hole, the 
time delay due to spacetime curvature (the Shapiro delay) 
is reduced. This means that for the most radially extended 
sources, the observed lag is less than that measured for the 
sources more confined to the region around the central black 
hole. 

Increasing the vertical extent of the source region (Fig. 
[5} again increases the range of possible ray paths from the 
source to the disc, giving the transfer function an extended 
tail and causing the measured lag to decay away more for the 
higher frequency components of the variability. Increasing 
the vertical extent of the source above a fixed bottom also 
moves the mean source location away from the accretion 
disc, increasing the mean lag time measured for the lowest 
frequency components. 



5 DILUTION OF SPECTRAL COMPONENTS 

The preceding calculations of reverberation time lags assume 
that the light curves of the directly observed continuum and 
reflected emission can be measured directly. In reality, how- 
ever, light curves are obtained in spectral bands which are 
merely dominated by either the continuum emission (the 
'hard' band at 1-4 keV) or the reflection from the accretion 
disc (the 'soft' band at 0.3-1 keV). This means that the mea- 
sured time lags will be 'diluted' as the measured 'reflection' 
will include promptly-arriving photons from the continuum 
and the 'continuum' will include late-arriving photons from 
the reflection. 

Fig. [6] shows the X-ray spectrum of 1H 0707-495 mod- 
elled as the sum of a direct power law component, reflection 
from the accretion disc and black body emission from the 
disc. Integrating the number of photons detected from each 
of the components over each energy band reveals that 0.3- 
1 keV band (dominated by the reflection component) con- 
tains a contribution from the primary continuum at 70 per 
cent of the number of photons detected in this band from the 
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Figure 6. The time-averaged X-ray spectrum of 1H 0707-495 ob- 
served using the EPIC pn detector on XMM Newton in October 
2008, modelled as the sum of a direct power law component, re- 
fl ection from an ionised accretion disc (using the REFLIONX model 
of lRoss fc Fabianll2005T) and black body emission from the disc. 



Table 2. Contributions from the spectral components of 1H 0707- 
495 from the January 2008 XMM Newton observation in the hard 
and soft bands from which reverberation time lags are calculated. 
Photon counts arc normalised to the component said to be dom- 
inating in that band. 



Component 


0.3-1 keV 


1-4 keV 


Power law continuum 


0.7 


1 


Reflection 


1 


0.6 


Black body 


0.4 


10" 4 



reflection. Likewise, the 1-4 keV band (dominated by the pri- 
mary continuum) contains a contribution from the reflected 
emission at 60 per cent of the number of photons detected 
in this band from the power law. We assume the flux in the 
thermal, black body emission to be constant and therefore 
not to contribute to the lag. 

To determine the effect of this dilution on the observed 
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Figure 5. |(a)| Transfer functions and |(b)| lag spectra for X-ray sources extended vertically from 2r g to varying heights above the 
plane of the accretion disc. In each case, the X-ray source is extended radially out to 35 r g . 
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Figure 7. The combination of the transfer functions for reflection 
from the accretion disc and direct emission from the X-ray source, 
in the ratio found in 1H 0707-495, to give the overall transfer 
function for the soft (0.3-1.0 keV) band. 



lag spectrum, a composite transfer function for the soft 're- 
flection' band was formed by summing the transfer functions 
due to reflection from the accretion disc and for the direct 
emission seen from the primary X-ray source, weighted such 
that the integrated photon counts in the constituents are 
split in the ratio observed in the spectrum (Fig. [?)). Like- 
wise, the transfer function in the hard 'primary' band is 
constructed as the sum of contributions from the direct 
continuum emission and reflection from the accretion disc. 
These combined transfer functions are used to obtain the 
light curves that would be observed in the hard and soft 
bands with Equations [6] and [7] Using the photon fluxes cal- 
culated above, the composite transfer functions are 



L source 



TfJ.3-l.0kcV = ^reflection + 0.7T SC 

T±.0-4.0kcV = + 0.6T rcrl 

cction 

Fig.[5]shows the effect of dilution of the spectral compo- 
nents in the 'primary' and 'reflection' energy bands on the 
lag spectrum. It is clear that when the contribution from 
the primary X-ray source is included in the 'reflection' band, 
the observed time lag is reduced by around 40 per cent com- 
pared to the 'actual' time lag that would be seen if it were 



Figure 8. The effect of dilution of spectral components in the 
'primary' and 'reflection' energy bands on the lag spectrum is 
illustrated by comparing the lag spectrum assuming a pure re- 
flection component, the lag spectrum obtained using a reflection 
transfer function combined with a transfer function of the primary 
X-ray source contributing to this energy band and a lag spectrum 
also accounting for the reflection component contributing to the 
'primary' energy band in the calculation of the observed arrival 
time of the direct emission. In each case, the source is taken to 
extend radially out to 35 r g and vertically between 2.0 and 2.5 r g . 



possible to directly observe the primary and reflected light 
curves. When, in addition to this, the contribution of the 
reflected component is included in the primary continuum 
band, the average time lag is further reduced. The observed 
time lag at 10 -3 Hz can be reduced by up to 75 per cent 
from the 'actual' time lag. 

It is evident that when interpreting the observed time 
lags between the primary continuum emission and reflection 
in terms of the light travel time to constrain the location of 
the X-ray source with respect to the reflector, it is important 
to interpret the observed 'hard' and 'soft' spectral bands as 
combinations of the primary and reflected emission, rather 
than simply using them as proxies for the respective com- 
ponents. 
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6 PROPAGATION EFFECTS 

When considering X-ray sources spanning radially out to 
around 35 r g , it is not physically meaningful to consider the 
source varying instantaneously through its entire extent, as 
has been assumed to this point where, in evaluating the 
transfer functions, all of the rays were taken to originate 
from the X-ray source at the same time. A fluctuation in 
the intensity of the primary X-ray source will propagate 
throughout its extent and in fact must do so in a time greater 
than the light travel time across the source region (which 
would be 35 r g /c for a source spanning 35 r g ). 

There are a number of possible scenarios for the propa- 
gation of fluctuations in the intensity of the primary X-ray 
source, which can arise due to changes in the density or 
energy/temperature of the particles in the corona. It is pos- 
sible that the fluctuations in the corona originate from the 
injection of energy by a process close to the central black 
hole and that fluctuations in the X-rays emitted from the 
corona originate near the centre and propagate outwards. 
Alternatively, it is possible that energy is injected into the 
corona from the surface of the accretion disc, either directly 
or through magnetic fields that are anchored to the accre- 
tion disc that accelerate the particles of the corona to high 
energies. In this case, one might expect that fluctuations 
arise from fluctuations in the accretion disc (for instance, 
variations in the mass accretion rate or in the magnetic flux 
densit y as field lines are 'acc re ted' inwards through t he disc 
as per lBeckwith et al.|[2009l or iMcKinnev et alj|20l3 ). Such 
fluctuations are likely to propagate inwards through the ac- 
cretion disc and will cause fluctuations in the corona that 
originate at the edge and move inwards. Finally, it could be 
that there is no ordered propagation of fluctuations through 
the corona, rather variations in the intensity of the primary 
X-ray source are stochastic in nature, occurring randomly 
throughout the source region (for example, they arise due 
to isolated magnetic flaring events with no change in the 
bulk properties of the corona over the time-scale of the ob- 
servations) . 

We investigate whether it is possible to discriminate be- 
tween these scenarios in the observed lag spectra by comput- 
ing transfer functions for the primary X-ray source and the 
reflection from the accretion disc in which the rays originate 
from the source at variable times. The influence of propaga- 
tion effects is investigated on the lag spectrum of an X-ray 
emitting region extending radially outward to 35 r g and be- 
tween 2 and 2.5 r g above the plane of the accretion disc and 
the full effects of dilution of the spectral components in the 
two energy bands are considered with the dilution fractions 
as found in 1H 0707-495. 

Random fluctuations in the X-ray source (if the changes 
in the observed X-ray luminosity originate in random flaring 
events throughout the source region whose net effect is to 
change the total count rate we observe; the variability arising 
probabilistically) are simulated by assigning a random start 
time to the rays up to a maximum value (here taken to be a 
multiple of the light travel time from the centre to the edge 
of the source region) which defines a characteristic timescale 
over which variations occur throughout the source region. 
These rays make up the transfer function, so in this model, 
constructed to give insight into how such a system would 
manifest itself in reverberation measurements, the random 
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Figure 9. The lag spectrum between the hard (1.0-4.0 keV, dom- 
inated by the primary continuum) and soft (0.3-1.0 keV, domi- 
nated by reflection from the acctrion disc) X-ray bands, taking 
into account the contributions of both spectral components to 
each band, allowing the rays from the primary X-ray source to 
have random start times up to a maximum (as a function of the 
radius of the source). The emitting region is located between 2.0 
and 2.5 r g above the accretion disc and extends radially out to 
35 rv. 



fluctuations are not the cause of the variability, rather given 
the observed variability, they dictate its spatial propagation 
through the source region and therefore how it is passed into 
the X-ray reflection. 

The lag spectrum arising from such a scenario, with in- 
creasing ranges in the ray start times in each transfer func- 
tion is shown in Fig. [9] Randomising the start times of rays 
through the source region results in the dilution of the mea- 
sured lag between the arrival of the primary and reflected 
X-rays, with the lag between the primary and reflected com- 
ponents tending to zero when the timescale of the random 
flaring events is much greater than the light travel time be- 
tween the X-ray source and the reflector. 

To simulate the effect of fluctuations in the source lumi- 
nosity propagating outwards from the centre of the source 
region, the start time of a ray is taken to be proportional to 
its originating radius within the region, where the constant 
of proportionality represents the speed at which fluctuations 
move radially through the source region, corresponding to 
dynamical timescales within the corona, and is the same 
for all rays. Likewise, the effect of fluctuations propagating 
inwards from the edge of the X-ray source is simulated by 
setting the start time of the rays proportional to their initial 
radial distance from the edge of the source region. 

The propagation of luminosity fluctuations inward from 
the edge of a defined source region delays the arrival time 
of the reflected radiation with respect to the primary con- 



tinuum. The resulting lag spectrum is shown in Fig. 10(a) 



As the speed at which fluctuations propagate through the 
source region is decreased, the measured lag increases and 
the lag at high frequencies is damped to zero over a broader 
frequency range. 

Conversely, Fig. 1 10(b) | shows theoretical lag spectra in 
which the fluctuation in source intensity propagates outward 
from the centre at varying speeds. While the propagation of 
fluctuations outward throughout the source region reduces 
the apparent time lag between the primary continuum and 
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reflection, the most notable effect is the turnover at low fre- 
quencies, < 10 -4 Hz, once the fluctuations propagate slower 
than approximately one tenth of the speed of light. In these 
situations, low frequency variations appear in the soft X- 
ray band, representing the 'reflection' before the hard band, 
representing the 'primary continuum.' 

Due to the propagation delay, X-rays from the outer 
parts of the source, which due to their larger area make up 
the majority of the flux received from the primary X-ray 
source, are delayed in their emission and thus their arrival 
at the observer. On the other hand, X-rays emitted from 
parts of the source closer to the black hole are more likely 
to be reflected from the accretion disc rather than observed 
directly in the continuum as rays passing closer to the black 
hole will be focussed on to the accretion disc by the strong 
gravity. This means that, comparatively, the reflected X-rays 
originate from parts of the source region closer to the black 
hole than the observed continuum and since the propagation 
time of the fluctuation is greater than the light travel time to 
the reflector, the average arrival time of the soft X-ray band 
at the observer is slightly earlier than the average arrival 
time of the hard band once the full contributions of the 
continuum and reflection in each band are accounted for 
(the dilution), resulting in the soft band slightly leading the 
hard band at the lowest variability frequencies (this effect 
is phase- wrapped out of the variability at higher frequencies 
where the soft band is seen to lag behind the hard band). 

If the speed at which the fluctuations propagate is fur- 
ther reduced, the low frequency variability in the observed 
hard X-ray band (dominated by the primary continuum) 
lags further behind that in the soft band while the high fre- 
quency part of the lag spectrum decays as the increasing 
lags cause these components to be phase-wrapped out. 



7 THE LAG SPECTRUM OF 1H 0707-495 

The lag spectrum obtai ned for the narrow line Seyfert 1 
galaxy, 1H 0707-495 (see IZoghbi et all |2010| . detailing the 
data reduction in the context of computing the lag spec- 
trum) between the 1.0-4.0 keV band (dominated by the 
power law continuum) 0.3-1.0 keV band (dominated by re- 
flection from the accretion disc) is compared to the theo- 
retical lag spectra computed from the reflection of an X-ray 
continuum arising from an extended source from the accre- 
tion disc. 

It is clear that the high frequency part of the lag spec- 
trum (> 10 -3 Hz) in which the soft 'reflection' band lags 
behind the hard 'primary' band by around 30 s is explained 
naturally by an X-ray source extending outward from the ro- 
tation axis to 35 r g and located around 2 r g above the plane 
of the accretion disc which, once the effects of dilution are 
taken into account matches the lag of 30 s at 10 -3 Hz and 
the decay of the lag towards higher frequencies (solid line 
in Fig. |8]) . This is consist ent with the extent of th e primary 
X-ray source derived by I Wilkins fc Fabianl |2012i ) from the 
observed emissivity profile of the reflected iron Ka line from 
the accretion disc, showing that the two independent tech- 
niques compliment one another in placing constraints on the 
properties of the source. 

It is assumed here that the X-ray source has a constant 
luminosity throughout its extent. In reality, it is likely that 



the X-ray emission will vary in luminosity based upon its 
location, with, perhaps, the emission being concentrated to- 
wards the centre of the source region. Assuming constant 
source luminosity, however, minimised the number of free 
parameters in this model and finding the X-ray source un- 
der this assumption has allowed the data itself to constrain 
the location of the bulk of the X-ray emission. This result is 
the vital first step in more detailed modelling of the struc- 
ture of the corona which must simultaneously satisfy the 
constraints from the observed lag spectrum, emissivity pro- 
file and energy spectrum. 

The simple model does not, however, explain the lag 
spectrum at low frequencies where the hard band is seen to 
lag behind the soft band. Such a turnover in the lag spec- 
trum may be, in part, accounted for by the propagation of 
luminosity fluctuations through the primary X-ray emitting 
region (in fact, fluctuations must propagate in a finite time 
through a spatially extended source). If fluctuations in the 
X-ray source luminosity propagate outward from the cen- 
tre of the source region at velocities of order one tenth of 
the speed of light, the low frequency variations in the soft 
band lead those in the hard band, while the high frequency 
variability shows the reverberation lag with the soft band 



lagging behind the hard, shown in Fig. 10(b) 



The hard X-ray band lagging behind the soft band at 
low frequencies may suggest a scenario in which energy is 
injected into the corona at the centre, close to the black 
hole, then propagates outwards. It should be noted that 
this model does not make any assumptions about the origin 
of the fluctuations in the X-ray source luminosity or their 
power spectral density (PSD). In fact, since the theoretical 
lag spectra have been computed by convolving the observed 
hard band light curve of 1H 0707-495 with the transfer func- 
tions corresponding to the X-ray source and the reflection, 
the observed PSD is accounted for inherently in the model. 
Neither are any assumptions made of the underlying energy 
transport mechanism between the accretion flow into the 
X-ray emitting region. Rather, we account for the observed 
lag spectrum of AGN in terms of the propagation of varia- 
tions at a characteristic velocity (analogous to a dynamical 
timescale) through an X-ray emitting region which is also 
constrained by independent, spectral measurements of the 
accretion disc emissivity profile, and the response of the re- 
flecting accretion disc. 

In January 2011, 1H 0707-495 was observed to fall into 
a low flux state in which the spectrum was almost en- 
tirely dominated by reflection from the accretion disc with 
almost no flux ob served from the power law continuum 
l|Fabian et al.ll201lh . The reflection spectrum from the ac- 
cretion disc was identified with a very steeply falling emis- 
sivity profile implying a small source concentrated close to 
the black hole. This can be reconciled with the energy that 
accelerates particles in the X-ray source originating from 
the central regions and then propagating outwards to form 
an extended source region. When the available energy de- 
creases, particles are still accelerated close to the black hole, 
though since less energy is available, it is not propagated 
outwards meaning the source is unable to extend, reducing 
the detected continuum emission as the emitting volume is 
smaller and the emission closer to the black hole tends to 
be focussed on to the accretion disc to be detected after 
reflection. 
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Figure 10. Lag spectra arising from fluctuations in the X-ray source luminosity propagating |(a)| inward from the edge of the 
defined source region and |(b)| outward from the centre. In each case, the lag spectra are shown for varying velocities at which 
the fluctuations propagate radially through the source region. The lag spectra are between the hard (1.0-4.0 keV, dominated by 
the primary continuum) and soft (0.3-1.0 keV, dominated by reflection from the acctrion disc) X-ray bands, taking into account 
the contributions of both spectral components to each band and the emitting region is located between 2.0 and 2.5 r g above the 
accretion disc and extends radially out to 35 r g . 



It is possible that the fluctuations in the luminosity 
of the primary X-ray source are due to fluctuations into 
the mass accretion rate onto the central black hole (vary- 
ing the energy that is released during the gravitational in- 
fall), with the fluctuations in accretion rate resulting in den- 
sity fluctuations spiralling inward s through th e accr etion 
disc. Such a mod e l is pr oposed bv iKotov et alJ (|200lh and 
lArevalo fc Uttlevl (|2006l ) in which the X-ray variability of 
Galactic black hole binaries is explained by fluctuations in 
the accretion disc propagating inwards, modulating the lu- 
minosity of the X-ray emitting region as they do so. This 
model does not account for how the accretion disc is cou- 
pled to the corona and only finds time lags in different energy 
bands of the emission from the primary X-ray source with 
measurements in energy bands both dominated by the con- 
tinuum. In the case of AGN, however, the dominant spectral 
component in the soft band is the reflection from the accre- 
tion disc with a characteristic lag spectrum resulting from 
reverberation of variability in the source. 

The simple reverberation from a fluctuation that propa- 
gates outwards through the X-ray source does not, however, 
reproduce the entirety of the hard lag which is observed to 
be longer. Other physical effects are likely to contribute to 
this. As well as fluctuations propagating through the source, 
the source properties are likely to vary throughout its extent, 
such as its density and energy per particle, meaning that the 
luminosity and spectrum (characterised by the power law 
slope) of X-rays emitted from different parts of the source 
may change. A greater Compton optical depth, up to unity, 
will result in the seed photons undergoing more scatterings 
and therefore reaching higher energies, producing a harder 
continuum spectrum. However, increasing the optical depth 
above unity will mean the photons undergo so many scat- 
terings, the net effect is that energy is distributed more be- 
tween the photons, resulting in a softer spectrum with more 
photons at lower energies. Likewise, increasing the average 
energy of the particles in the corona will produce a harder 



spectrum where more ene rgy is transferred to the photons 
(see, e.g., lTitarchukllT99i ). 

The observed hard lag at low frequencies can result from 
the X-ray source itself if the fluctuations reach parts of the 
source producing a harder spectrum after they reach parts 
of the source producing a soft spectrum. This will mean at 
low frequencies, variability in the hard X-ray band will lag 
behind that in the soft band. If this effect is to add to the 
hard lag produced intrinsically by the propagation of fluctu- 
ations through the source, the central regions of the X-ray 
source produce a softer power law continuum spectrum than 
the outer parts. This is possible if the average energy of the 
particles in the corona is greater in the outer parts. If it 
is only the higher energy particles that are able to travel 
out further from the central part into which the energy is 
injected from the accretion flow with the lower energy parti- 
cles remaining in the centre, the average energy of particles 
in the outer parts of the source region will be greater, pro- 
ducing a harder spectrum here. Alternatively, a softer power 
law continuum could be produced by the innermost parts of 
the source if the density is much greater here such that the 
optical depth is greater than unity (such that the increasing 
density produces a softer spectrum), though it is unlikely a 
Compton-thick source extends further out than the inner- 
most stable orbit of the accretion disc as this corona would 
obscure the reflected emission seen from the disc. 

Finally, it is possible that the fluctuations in source lu- 
minosity propagate inwards from the edge of the source re- 
gion and the inner parts of the X-ray source produce a harder 
power law continuum as the Compton depth and particle en- 
ergy increases towards the centre, however for this to be the 
case, the variation in spectral slope but be sufficient across 
the source that the intrinsic time lag between the energy 
bands within the source is able to counteract the increased 
soft lag (the soft, reflection band lagging behind the hard, 
continuum band) that results from the fluctuations propa- 
gating inwards. 

The applicability to 1H 0707-495 of a simple model con- 
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sisting of two power law continua with different slopes with 
the hard power law laggi ng behind the soft reproducing the 
hard lag is considered bv lKara et al.l (|2012h . Detailed mod- 
elling of the lag spectra arising from such a scenario would 
require detailed analysis of the spectral contributions that 
could come from power law continua with different slopes 
and will be considered more fully in future work. 



8 A SECOND EXAMPLE: IRAS 13224-3809 



iFabian et aD l|2012t ) obtained the lag spectrum of the nar- 
row line Seyfert 1 galaxy IRAS 13224-3809 from a long ob- 
servation of the X-ray emission from this source with XMM 
Newton. The lag spectrum is of the same form is that of 
1H 0707-495 with a reverberation lag (the soft band corre- 
sponding to reflection from the accretion disc lagging behind 
the hard band corresponding to the directly observed con- 
tinuum) of around 100 s at a frequency of around 4 x 10 -3 Hz 
and the turnover at low frequencies with the hard band lag- 
ging behind the soft. 

The longer observed lag (once dilution has been taken 
into account) and the shift of the form of the lag spectrum 
to lower frequencies is consistent with the greater mass of 
the central black hole in IRAS 13224-3809 than in 1H 0707- 
495 (around 10 7 M©) and with the X-ray source being more 
compact. The emissivity profile of the accretion disc implies 
the X-ray source is contained within the innermost 10 r g , 



which, from Fig. 4(b) implies in itself a shift of the rever- 
beration lag to lower frequencies and a longer reverberation 
lag than for the more extended source in 1H 0707-495 since 
there are fewer possible ray paths from an extended source 
and the rays pass closer to the central black hole, increasing 
the influence of the Shapiro delay. 



9 CONCLUSIONS 

In order to interpret observed lag spectra and reverbera- 
tion lags between the primary continuum and reflected X-ray 
emission in AGN, theoretical lag spectra were calculated for 
reflection of the X-ray continuum from the accretion disc 
by computing transfer functions in general relativistic ray 
tracing simulations which were convolved with the observed 
hard band X-ray light curve of 1H 0707-495. Lag spectra 
were computed for isotropic point sources at varying heights 
as well as a variety of spatially extended sources. 

The Shapiro delay of light travelling close to the cen- 
tral black hole as a significant effect on the travel time of 
rays from the source to the reflector and to the observer. 
Therefore, when inferring the location and properties of the 
primary X-ray source from measurements of reverberation 
lags, it is important to account for general relativistic effects 
rather than directly equating the measured lag to a light 
travel time and the equivalent distance in flat, Euclidean 
space. 

When interpreting observed reverberation lags, it is nec- 
essary to consider the contribution of all spectral compo- 
nents to the 'hard' and 'soft' bands rather than taking them 
to be direct proxies for the primary and reflected radiation. 
The 'hard' and 'soft' X-ray bands are seen to have signif- 
icant contributions from both the primary continuum and 



reflected spectral components in AGN such as 1H 0707-495 
which causes the observed time lag between the arrival of the 
primary continuum and reflection components to be diluted, 
reducing the observed lag by up to 75 per cent compared to 
the 'intrinsic' light travel time delay. It is therefore vital to 
have adequate spectral models to understand the contribu- 
tions to the energy bands when interpreting reverberation 
lags. 

By comparing the observed lag spectrum of 1H0707- 
495, we are able to constrain the primary X-ray continuum 
to originating from a source region extending radially out- 
wards to around 35 r s from the rotation axis, located around 
2r g above the plane of the accretion disc. This is consistent 
with constraints on the location and extent of the X-ray 
source derived independently from the observed X-ray re- 
flection emissivity profile of the accretion disc. 

The propagation of luminosity fluctuations through an 
extended X-ray source is considered, since it is unphysical to 
think of the luminosity of a spatially extended source varying 
instantaneously. It is demonstrated that if the fluctuations 
in luminosity are allowed to propagate radially outward from 
the centre of the source at a velocity (corresponding to a dy- 
namical timescale in the corona) of order a tenth of the speed 
of light, the resulting lag spectrum turns over at low frequen- 
cies, < 10 -4 Hz and variability in the soft band (dominated 
by reflection) leads that in the hard band (dominated by the 
primary continuum) as is seen in AGN such as 1H 0707-495. 
This suggests that the hard lag could arise, in part, from 
a scenario in which energy is injected into the corona from 
the accretion flow very close to the central black hole and 
then propagates outwards into the spatially extended source 
region. 

In this work, no assumptions are made as to the under- 
lying physical processes through which energy is released 
from the accretion flow into the X-ray emitting corona, 
rather constraints are placed upon the properties the X- 
ray source must have to be consistent with observations in 
the context of the reverberation off an accretion disc sur- 
rounding a black hole in general relativity. These findings 
may be compared to more physical models of coronae or in- 
deed rays may be traced in the same framework from X-ray 
sources produced by these models to constrain the under- 
lying physical processes. It is argued that the energetics of 
the X-ray emitting corona must be dominated by magnetic 
fields and corona composed of magnetic flux loops emerg- 
ing from the accretion disc due tot he magneto-rotational 
instability have been studied extensively by, among others, 
lUzdenskv fc Goodman! f2008). They describe the corona as 
a statistical ensemble of magnetic flux loops in (magnetic) 
pressure equilibrium with their surroundings, allowing them 
to undergo reconnection and predict that the magnetic en- 
ergy density should be concentrated at low heights above 
the accretion disc (within a few times the scale height of 
the disc). Assuming a thin accretion disc with aspect ratio 
H/r ~ 0.1, this is roughly consistent with our findings from 
the reverberation lags that the bulk of the X-ray source is 
located within a few gravitational radii above the accretion 
disc, though more detailed theoretical work tailored to these 
observations would be required to place any firm constraints 
on the underlying process powering the corona. 
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